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Abstract 

Recently, researchers have apphed genetic algorithms (GAs) to address some problems 
in quantum computation. Also, there has been some works in the designing of genetic 
algorithms based on quantum theoretical concepts and techniques. The so called Quan- 
tum Evolutionary Programming has two major sub-areas: Quantum Inspired Genetic 
Algorithms (QIGAs) and Quantum Genetic Algorithms (QGAs). The former adopts 
qubit chromosomes as representations and employs quantum gates for the search of 
the best solution. The later tries to solve a key question in this field: what GAs will 
look like as an implementation on quantum hardware? As we shall see, there is not a 
complete answer for this question. An important point for QGAs is to build a quantum 
algorithm that takes advantage of both the GA and quantum computing parallelism 
as well as true randomness provided by quantum computers. In the first part of this 
paper we present a survey of the main works in GAs plus quantum computing including 
also our works in this area. Henceforth, we review some basic concepts in quantum 
computation and GAs and emphasize their inherent parallelism. Next, we review the 
application of GAs for learning quantum operators and circuit design. Then, quantum 
evolutionary programming is considered. Finally, we present our current research in 
this field and some perspectives. 
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1 Introduction 

Our aim in this paper is two-fold. Firstly, we review the main works in the application of 
Genetic Algorithms (GAs) for quantum computing as well as in the Quantum Evolutionary 
Programming. Secondly, based on this review, we offer new perspectives in the area which 
are part of our current research in this field. 

In the last two decades we observed a growing interest in Quantum Computation and 
Quantum Information due to the possibility to efficiently solve hard problems for conven- 
tional computer science paradigms. Quantum computation and quantum information en- 
compass processing and transmission of data stored in quantum states (see POj and refer- 
ences therein). In these fields, the computation is viewed as effected by the evolution of a 
physical system, which is governed by unitary operators, according to the Laws of Quan- 
tum Mechanics |19j . The basic unity information is the qubit, the counterpart in quantum 
computing to the classical 0—1 bit. Quantum Computation and Quantum Information 
explore quantum effects, like quantum parallelism, superposition of states and entanglement 
in order to achieve a computational theory more efficient than the classical ones. This has 
been demonstrated through quantum factoring and Grover's algorithm for database search 
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On the other hand, Genetic Algorithms (GAs) is a rapidly expanding area of current 
research. They were invented by John Holland in the 1960s Simply stated, GAs are 
stochastic search algorithms based on the mechanics of natural selection and natural genetics 
Pl IlKL [T5] . They have attracted people from a wide variety of disciplines, mainly due to 
its capabilities for searching large and non-linear spaces where traditional methods are not 
efficient 

From this scenario, emerge the application of genetic algorithms for quantum compu- 
tation as well as evolutionary programming based on quantum theoretical concepts and 
techniques. Despite of the fact that there are few works in these subjects yet, it is an 
exciting area of research in the field of evolutionary computation. 

When applying GAs, people are attracted by their capabilities for searching a solution 
in non-usual spaces. That is way people investigate the application of GAs for learning 
quantum operators |S1|S| and in the designing of quantum circuits [5111^0^) . 

Those works rely on the fundamental result for quantum computing that all the compu- 
tation can be expanded in a circuit which nodes are the universal gates |19j . These gates 
offer an expansion of an unitary operator U that evolves the system in order to perform 
some computation I14j . Thus, we are naturally in the face of two classes of problems: 
(1) Given a set of functional points S = {{x, y)} find the operator U such that y = U ■ x; (2) 
Given a problem, find a quantum circuit that solves it. The former was formulated in the 
context of GAs for learning algorithms [HI El while the latter through evolutionary strategies 
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In we proposed a method based on genetic algorithms to learn linear operators. The 
method was applied for learning quantum (unitary) operators. It was demonstrated that it 
overcomes the limitations of the work proposed by Dan Ventura j28J, which resembles basic 
methods in neural networks. 

For the second class of problems, we found three schemes outlined by Spector 1^ , based 
on the traditional tree-based genetic programming stackless and stack-based linear 
genome. Moreover, in another scheme proposed by Williams and Gray |29j . an unitary 
matrix for a known quantum circuit is used to find possible alternative circuits. 

More close to our actual research are the works by Rubinstein and Yabuki J31 . The 
production of entangled states and the quantum teleportation were the target problems in 
these works. In each gate is encoded through a structure (gate structures) and a quantum 
circuit is represented as a list of such structures. Then, genetic operators (crossover and 
mutation) are applied in order to evolve a randomly chosen initial population of circuits. 
A special fitness function was also proposed and the scheme applied to find circuits for 
entangled states production. Yabuki and Iba follow a similar philosophy in but changing 
the circuit representation. The quantum teleportation was the focused problem. In this case, 
authors performed a circuit optimization; that is, they start with a seed circuit that had 
eleven gates and obtained another one, with just eight gates, but that performs the same 
computation. That is a variant of the second class of problems. 

By 1996, Narayan and Moore |17| introduced a novel evolutionary computing method 
where concepts and principles of quantum computing are used to inspire evolutionary strate- 
gies. It was the first attempt towards quantum evolutionary programming. The basic ap- 
proach is inspired on the multiple universes view of quantum theory: each universe contains 
its own population of chromosomes. The populations in each universe obey identical classi- 
cal rules and evolve in parallel. However, just after classical crossover within each universe, 
the universes can interfere with one another which induces some kind of crossover involving 
the chromosomes. 

Narayan and Moore's method depends on non-standard interpretations of Quantum Me- 
chanics. The lack of a more formal analysis of the physical concepts used brings difficulties 
to make the correlation between the physics and the genetic algorithm itself. Consequently 
it does not offer clues of the advantages of a quantum implementation for GAs within the 



current implementations of quantum computers j2U| . Thus, we are not going to consider it 
in the following sections. 

In jllj we found another Quantum Inspired Genetic Algorithm (QIGA) which relies 
on usual methods of Quantum Mechanics. It is characterized by principles of quantum 
computing, including concepts of qubits and superposition of states, as well as quantum 
operators to improve convergence. An important consequence of this work is to emphasize 
that the application of quantum computing concepts to evolutionary programming is a 
promising field. 

Rylander et al. kept this idea and sketched out a Quantum Genetic Algorithm (QGA) 
which takes advantage of both the quantum computing and GAs parallelism. Despite of the 
lack of a mathematical explanation about the physical realization of the algorithm, we will 
show that its philosophy is close to an implementation in current experimental architectures 
for quantum computers j20| . The key idea is to explore the quantum effects of superposition 
and entanglement to create a physical state that store individuals and their fitness. When 
measure the fitness, the system collapses to a superposition of states that have that observed 
fitness. Starting from this idea, we propose in section [S] a QGA which can take advantage 
of both quantum computing and GAs paradigms. We present its physical foundations and 
discuss its advantages over classical GAs. 

This paper is organized as follows. Firstly, we give the necessary background in quantum 
computation fsection |2.1|) and genetic algorithms fsection |2.2|I . Then, the review sections 
start. We begin with the GAs applications in quantum computation. So, in section [3. II we 
describe our work on GAs for learning quantum operators. Section lX^ presents the works on 
GAs for circuit design. Following, in sectional we offer the review of quantum evolutionary 
approaches. The QIGA proposed in is presented on section im We end the review by 
presenting on section the QGA proposed in |2S]. In sectional we discuss the considered 
methods and describe some issues. In particular, we propose a physical model for the QGA. 
Finally, we present the conclusions on section IHl 

2 Background in Quantum Computation and GAs 
2.1 Quantum Computation 

In practice, the most useful model for quantum computation is the Quantum Circuit one 
\W\ 12] • The basic information unit in this model is a qubit , which can be considered 
a superposition of two independent states | 0) and | 1) , denoted by | t/;) = ao I 0) + ai | 1), 
where aoiQ^i are complex numbers such that \ao\^ + |q;i|^ = 1. They are interpreted as 
probability amplitudes of the states | 0) and | 1). 

A composed system with n qubits is described using TV = 2" independent states obtained 
through the tensor product of the Hilbert Space associated with each qubit. Its physical 
realization is called a quantum register. The resulting space has a natural basis that can be 
denoted by: 



where we are using the Dirac notation for vectors in Hilbert spaces. 

This set can be indexed by|i); i — 0,1, N — 1. Following the Quantum Mechanics 
Postulates , the state | i/') j of a system, in any time t, can be expanded as a superposition 
of the basis states: 
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Entanglement is another important concept for quantum computation with no classical 
counterpart. To understand it, a simple example is worthwhile. 

Let us suppose that we have a composed system with two qubits. According to the above 
explanation, the resulting Hilbert Space has N — 2"^ independent states. 

Let the Hilbert Space associated with the first qubit (indexed by 1) denoted by Hi 
and the Hilbert Space associated with the second qubit (indexed by 2) denoted by H2. The 
computational basis for these spaces are given by: {| 0)i, | l)i} and {| 0)2, | 1)2}, respectively. 
If qubit 1 is in the state | -0)1 = aio I 0)i + an | l)i and qubit 2 in the state | -0)2 = 020 | 
0)2 + 021 I 1)2, then the composed system is in the state: | ^) =\ V')i<8i | ■0)2, explicitly given 
by: 

I V) = I «)l«) I j)2- (3) 
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Every state that can be represented by a tensor product | \ 0)2 belongs to the 

tensor product space Hi (g) H2. However, there are some states in Hi (g) H2 that can not be 
represented in the form | 'ip)i'^ \ V')2- They are called entangled states. The Bell state (or 
EPR pair), denoted by | /3oo)j is a very known example: 

I /3oo) = ^ (I 0)1® I 0)2+ I 1)1® I 1)2). (4) 

Trying to represent this state as a tensor product | 0)ig) | 0^)2, with \ -1^)1 £ Hi and 
I 0)2 G H2, produces an inconsistent linear system without solution. 

Entangled states are fundamental for teleportation [71^]. In recent years, there has been 
tremendous efforts trying to better understand the properties of entanglement, not only as 
a fundamental resource for the Nature, but also for quantum computation and quantum 
information 

The computation unit in quantum circuits's model consists of quantum gates which are 
unitary operators that evolve an initial state performing the necessary computation. A 
quantum computing algorithm can be summarized in three steps: (1) Prepare the initial 
state; (2) A sequence of (universal) quantum gates to evolve the system; (3) Quantum 
measurements . 

From quantum mechanics theory, the last stage performs a collapse and only what we 
know in advance is the probability distribution associated to the measurement operation. 
So, it is possible that the result obtained by measuring the system should be post-processed 
to achieve the target (quantum factoring (Chapter 6 of 21 ) is a nice example). 

More formally, the measurement in quantum mechanics is governed by the following 
postulate |19j . Quantum Measurements are described by a collection {Mm} of measurement 
operators satisfying {J2m ^'^m^'^m = 1) acting on the state space of the system being mea- 
sured. If the state on the system is | t/i), given by expression immediately before the 
measurement then the probability that result m occurs is given by: 

p{m)^<^\MlM^\^), (5) 
and the state of the system just after the measurement is: 

M„ 10) 

I V) after = , (6) 
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The expression (jSJ is the mathematical description of the collapse due to the measure. 
A simple but important example is the measurement of the state | ip) , given by expression 
((21) in the computational basis. This is the measurement of the system with N outcomes 



defined by the operators {| 0) < |, | 1) < 1 |, | 2) < 2 |, | - 1) < iV - 1 |} . Thus, ap- 
plying equations (ISl)-® we find: 

p{m) = \amf, I V')Q/ter = -p^ I m = 0, . . . , TV - 1 . (7) 

I 

Thus, the collapse becomes more evident. 

Quantum parallelism is another fundamental feature of quantum computing. To better 
explain it, let us take the Hadamard Operator (H) defined by: 

H|0) = ^, (8) 

Now, we shall present another quantum operator which will be useful in the following 
sections. Suppose / : {0, 1} — > {0, 1} a binary function with a one-bit domain. Now, define 
a quantum operator U f : Hi ® H2 ~* Hi ® H2 such that: 

Uf\x)®\y)^\x)®\y®fix)), (9) 

where the symbol © means addition modulo 2. If we take the state | -0) =| 0)® | 0) and 
apply the Hadamard operator over the second qubit (/ (g) H), followed by Uf we obtain: 



il^HU 0>« I 0) = U, I 0) « =1 0) « ([imiml^ ^ ,10, 

This is a remarkable result because it contains information about / (0) and / (1) . It is 
almost as if we have evaluated / (x) for two values of x simultaneously! This feature is 
known as quantum parallelism and can be generalized to functions on an arbitrary number 
of bits. 

To simplify notations, we will represent the tensor product | i)<8' | j) by | ij) or | z) | j), 
in what follows. 

A fundamental result for quantum computing is that any unitary matrix U which acts 
on a d— dimensional Hilbert space can be represented by a finite set of unitary matrices 
( Universal Gates) which act non-trivially only on lower subspaces. The Hadamard operator 
defined in expression JHJ is an universal quantum gate. Another one is the CNOT operator, 
defined as follows: 



CNOT I 00)=|00); C7VOT | 01) =| 01); 
CNOT I 10) =1 11); CNOT \ 11) =| 10); 



(11) 



that is; the action of the operator is such that if the first qubit {control qubit) is set to zero, 
then the second qubit {target) is left unchanged. Otherwise, the target qubit is flipped. 
Other universal quantum gates can be found in |19| . 

Given an operator in a Hilbert space, we can take its action over the computational basis 
to get a matrix representation. For the Hadamard and CNOT gates, we have the following 
representations: 
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10 
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according to expressions ^ and lfTT|l , respectively. 



2.2 Evolutionary Computation and GAs 

In the 1950s and the 1960s several computer scientists independently studied evolutionary 
systems with the idea that evolution could be used as an optimization tool for engineering 
problems. The idea in all these systems was to evolve a population of candidate solutions 
for a given problem, using operators inspired by natural genetic and natural selection. 

Since then, three main areas evolved: evolution strategies, evolutionary programming, 
and genetic algorithms. Nowadays, they form the backbone of the field of evolutionary 
computation [TCllTj. 

Genetic Algorithms (GAs) were invented by John Holland in the 1960s TS'. Holland's 
original goal was to formally study the phenomenon of adaptation as it occurs in nature 
and to develop ways in which the mechanism of natural adaptation might be imported into 
computer systems. In Holland's work, GAs are presented as an abstraction of biological 
evolution and a theoretical framework for adaptation under the GA is given. Holland's GA 
is a method for moving from one population of chromosomes to a new one by using a kind 
of natural selection together with the genetic-inspired operators of crossover and mutation. 
Each chromosome consists of genes (bits in computer representation), each gene being an 
instance of a particular allele {O or 1). 

Traditionally, these crossover and mutations are implemented as follows ^] . 

Crossover: Two parent chromosomes are taken to produce two child chromosomes. Both 
parent chromosomes are split into left and a right subchromosomes. The split position 
(crossover point) is the same for both parents. Then each child gets the left subchromosome 
of one parent and the right subchromosome of the other parent. For example, if the parent 
chromosomes are Oil 10010 and 100 11110 and the crossover point is between bits 3 and 4 
(where bits are numbered from left to right starting at 1), then the children are 01111110 
and 100 10010. 

Mutation: : When a chromosome is taken for mutation, some genes are randomly chosen 
to be modified. The corresponding bits are flipped from to 1 or from 1 to 0. 

These operations reveal the fact that GAs are inherently parallel algorithms. GAs work 
by discovering the most adapted chromosomes, emphasizing, and recombining their good 
"building blocks" through operations that can be easily performed in parallel. This has 
been explored in many works over the GA literature |31 112j . 

Next we present a generalization of the 0—1 case, in which the alleles are real parameters 
It belongs to the class of real-coded Genetic Optimization Algorithms and will be used 
in section im 

Genetic Optimization Algorithms are stochastic search algorithms which are used to 
search large, non-linear spaces where expert knowledge is lacking or difficult to encode and 
where traditional optimization techniques fall short [S]. 

To design a standard genetic optimization algorithm, the following elements are needed: 

(1) A method for choosing the initial population; 

(2) A "scaling" function that converts the objective function into a non-negative fitness 
function; 

(3) A selection function that computes the "target sampling rate" for each individual. 
The target sampling rate of an individual is the desired expected number of children for 
that individual. 

(4) A sampling algorithm that uses the target sampling rates to choose which individuals 
are allowed to reproduce. 

(5) Reproduction operators that produce new individuals from old ones. 

(6) A method for choosing the sequence in which reproduction operators will be applied 
For instance, in |3()j each population member is represented by a chromosome which is the 

parameter vector x — [xi, X2, Xm) G Si™, each component Xi being a gene. Consequently, 
alleles are allowed to be real parameters. Thus, some care should be taken to define these 



operators. 

Mutations can be implemented as a perturbation of the chromosome. In [30] , the authors 
chosen to make mutations only in coordinate directions instead of to make in 3?™ due to 
the difhcult to perform global mutations compatible with the schemata theorem (it is a 
fundamental result for GAs |13l l^). 

Besides, the crossover in 3?™ may also have problems. Figure Q illustrates the difficult. 
The ellipses in the figure represent contour lines of an objective function. A local minimum 
is at the center of the inner ellipse. Points {xi,yi) and (x2,?/2) are both relatively good 
points in that their function value is not too much above the local minimum. However, if we 
implement a traditional- like crossover (section 12.21) we may get points that are worse than 
their parents. 




Figure 1: Crossover can generate points out of the attraction region. 

To address this problem, in '30' was proposed another form of reproduction operator 
that was called linear crossover. From the two parent points pi and p2 three new points are 
generated, namely: 

1 , ,31 13 

^{pi+P2); -P1--P2; ~-pi + -p2. 

The best two of these three points are selected. 

Inspired on the above analysis we propose the algorithm of section 13.11 to learn a linear 
operator from a set S of example functional points. 



3 Applying GAs for Quantum Computing 

3.1 GA for Learning Operators 

Let us suppose that we do not know an operator F : V V but, instead, we have 
a set of functional points S = {(| Xi), I ^i)) P \ Xi) =\ * = 0, 1, n — 1} , where 

dim iy) = n, also called the learning sequence. We can hypothesize a function G such that 
\\G I Xi)~ I V'i)!! — (as usual, ||| v)\\ — \/ {v \ v) is the norm induced by the inner product). 

In this section, we present our general learning algorithm, based on GAs, to find G 0. 
This work was motivated by Dan Ventura's algorithm for learning quantum operators |28j . 
which resembles basic methods in neural networks 0. Our GA method has a range of 
applications larger than that one of Dan Ventura's learning algorithm. This is the main 
contribution of the work described next. 

Following [30], each population member (chromosome) is a matrix A e 3?"^", and alleles 
are real parameters (matrix entries). The alleles are restricted to [—1,1], but more general 
situations can be implemented. 



The initial population is randomly generated. Once a population is obtained, a fitness 
value is calculated for each member. The fitness function is defined by: 

fitness {A) ^exp {-error (A)); A e 3?"^", (13) 

where the error function is defined as follows. Let the learning sequence S = 
{(I Xi)J V"*)) ; i = 1, ...,to} , then: 

error (A) ^ ^ P | ^^)_ | ^^)|| (14) 

n ■ m ^ — ' 

i=l 

where ||a;||-^ denotes the 1-norm of a a; = (xi, defined by: ||a;||-^ — \xi\ + ... + \xn\ . 

Once the fitness is calculated for each member, the population is sorted into ascending 
order of the fitness values. Then, the GA loop starts. Before enter the loop description, 
some parameters must be specified. 

Population Size: Number of individuals in each generation {N). 

Elitism: It might be convenient just to retain some number of the best individuals of 
each population (members with best fitness) {Ne). The other ones will be generated through 
mutation and/or crossover. This kind of selection method was first introduced by Kenneth 
De Jong |16| and can improve the GA performance. 

Selection Pressure: The degree to which highly fit individuals are allowed many off- 
springs (Ps). For instance, for a selection pressure of 0.6 and a population with size N, 
we will get only the 0.6 * N best chromosomes to apply genetic operators. 

Mutation Number: Maximum number of alleles that can undergo mutation (Nm). 
We do not choose to make mutations (implemented as perturbations) in 3?"^", likewise in 
j3()| . Instead, we randomly choose some matrix entries to be perturbed. 

Termination Condition: Maximum number of generations (Ngen). 

Mutation and Crossover Probabilities: Pm and Pc, respectively. 

The crossover is defined as follows. Given two parents A = [aij] and B ~ [bij] the 
following steps are performed until two offspring Ci, C2 are generated: (1) Randomly choose 
one of the parents; (2) Take a (ij) matrix entry and puts its value on Cij. Go to step (1). 

The mutation is implemented as a perturbation of the alleles. Thus, given a member A, 
the mutation operator works as follows: A —^ A + A; where A is a perturbation matrix. 
The mutation number establishes the quantity of non-null entries for A. They are defined 
according to the mutation probability and a pre-defined Perturbation Size, that is, a 
range [a, b] € 3?, such that a < Ay < b. 

Once the above parameters are pre-defined and the input set (5) is given, the GA al- 
gorithm proceeds. In the following pseudo-code block, P {t) represents the population at 
the interaction time t and N is its size. Ngen is the maximum number of generations al- 
lowed, the procedure Evaluate_Sort calculates the fitness of each individual and sort the 
chromosomes into ascending order of the fitness values. The integer A^e defines de elite 
members, the parameter Ps G [0, 1] defines the selection pressure and Nm the number of 
matrix entries that may undergo mutations. 

Procedure Learning-GA 
t ^ 0; 

initialize P (t); 
while(t < Ngen) do 

EvaluateJSort {P [t — 1)); 

Store in P {t) the A^e best members of P (t — 1); 
Complete P {t) by crossover and mutation; 
end while 



3.1.1 Experimental Results 



Firstly, we analyze the behavior of the GA learning algorithm for the same example presented 
in [2H|- The set S is given by: 
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and the target is the Hadamard Transform H, defined by expression JHJ. 

The GA result over 25 runs was always the correct one. The set of parameters is given 
in the first line of Tabled The perturbation size is given by [0.001, 0.1]. 

Matrix Ngen N Pc Pm Ps Ne Nm 
2 x 2 100 200 0.85 0.95 0.30 30 1 
2 x 2 200 200 0.85 0.95 0.30 30 1 



Table 1: Parameters for the examples of this section (definitions on section m|l . 

Figure El shows the error evolution over the 25 runs. We collect the best population 
member (smallest error) for each run and take the mean value, for each generation, over the 
25 runs. It suggests that the algorithm gets closer the solution fast but takes much more 
time to achieve the target. Indeed, this behavior was observed for all experiments reported 
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Figure 2: Error evolution, over 25 runs, for example given by expression {T^J. 
Dan Ventura's algorithm gives also the correct result for this example, as reported in 
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Now, let us take the following set S: 
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In this set, the input vectors are not orthonormal ones. Thus, as we demonstrate in [H], 
if we apply Ventura's algorithm we get a result which is far from the target (see [S| for more 
details). But, our GA algorithm was able to deal with this case. The operator to be learned 
is the Hadamard Transform, as before. 

The second line of Table shows the parameters used. We would like to keep all pa- 
rameters unchanged but the number of generations (Ngen) had to be increased to achieve 
the correct result. This point out that our GA method may be sensitive to the fact that 
the set {| xo)i I Xi)} i^ot orthonormal basis, despite that it learns correctly. Additional 
examples must be performed in order to verify this observation. 

The mean error evolution shows a behavior which is similar to the first example. It 
decays fast but takes some time to become null. 

Table ^ shows that we could keep the crossover and mutation probabilities unchanged 
for both these experiments (more results are presented on It is desired because it may 
indicates some parameter stability. Moreover, the clock time for one run is very acceptable 
(< 0.04 seconds). 



3.2 GA for Quantum Circuit Design 

Despite of the scientific and technological importance of quantum computation, few quantum 
algorithms faster then the classical ones have been discovered. Shor's algorithm for quantum 
factoring. Grove's quantum search and Deutsch-Jozsa algorithm are basically the known ones 

oiini. 

This is due to the fact that the generation of such algorithms or circuits is difhcult for 
a human researcher. They are unintuitive, mainly due to quantum mechanics features like 
entanglement and collapse fsection l2.1|l . That is way researchers have investigated the use 
of stochastic search techniques, such as genetic programming and genetic algorithms to help 
in this task. 

Among the main works in this subject |31[ I23L I29L I27| . those ones proposed in j2I] and 
PS) are closer to our current research in this field. 

The paper presents a new representation and corresponding set of genetic operators 
for a scheme to evolve quantum circuits. A quantum circuit is represented as a list of 
structures (gate structures), where the size of a circuit (number of gates) can vary up to a 
pre-defined maximum number. Each gate structure contains the gate type, which is one in an 
allowable set of gate types including the usual Identity, CNOT, Hadamard and measurement 
operators, and a binary string for each category of qubit and parameter for that gate (see 
Figure El . 



(a) 
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Figure 3: (a) Quantum Circuit. From the left to right: Hadamard {H), CNOT and the 
Rotation Gate U (9). (b) The corresponding representation. 



The crossover and mutation operators are defined as follows. Crossover operates on all 
the levels of an individual's structure: the gates, qubit operands and parameter type. Gate 



crossover between two parent circuits consists of picking a gate from each parent at random, 
and then swapping aU gates between the parents after these two points. Crossover between 
binary strings representing parameters can only occur between strings of hke category, and 
proceeds in the same way as the crossover operator for the fixed length GA: pick a crossing 
point, and then swap bit values between the two strings after the point. 

Mutation happens in the gate level. A gate is mutated by replacing it with a new one 
randomly selected. Mutation is performed with small probability (typically 0.001) because in 
|23 l mutation is considered more an insurance against loss of important building blocks than 
a fundamental search procedure. Such viewpoint may be changed for circuit optimization. 
We shall return to this point ahead. 

An error function is defined to compare the stored state(s) with the desired one(s): 
Given a set of cases consisting of input states and desired outputs, the error is defined by 
error = J2i j Wij ~ where we take for each case i, the sum of the magnitudes of the 
differences between the probability amplitudes of the desired result di and that obtained one 
(Tj. A fitness function is constructed based on this error The production of entangled 
states was the focused application. 

Reference '31' is another proposal in the application of GAs for quantum circuit design. 
The philosophy is similar to the one presented above. 

The case study is the teleportation circuit 0. Quantum teleportation is a technique 
by which a quantum state can be transported from one point to another through non-local 
interactions To illustrate the steps involved in quantum teleportation let us consider 
that there are two friends, Alice and Bob. 

Imagine that Alice wants to deliver a qubit | /) = p | 0) + g | 1) to Bob, who lives far 
apart. Alice does not know p and q. Moreover, she can not have access to these values by 
measuring the qubit because, according to the quantum mechanics postulates, the system 
will collapse to a state p/ \p\ \ 0) or q/ \q\ \ 1). Once measurement is an irreversible task, 
the information would be lost. 

The scheme to solve this problem comprises the quantum teleportation. Its physical 
basis was proposed by Bennet et al. 15, followed by Brassard [5], who proposed the quantum 
circuit of Figure 01 for teleporting a single qubit. 
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Figure 4: Brassard's circuit for teleporting a single qubit. 

In the circuit of FigureQ] we have the following components: quantum gates L, i?, S", T, 
given by: 
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(18) 



measurement operator M and the CNOT gate, already defined in expression (|ll|l . and 
represented like in Figure 13 



The circuit has 3 qubits, namely qubit 0, 1 and 2, from the bottom to the top of the 
circuit. The method explores the concept of entanglement by using the EPR state given 
in equation The first part of the circuit is used to create that kind of entangled state. 
In this operation, the zeroth and first qubits are affected. The input state is given by the 
tensor product | /) | 0) | 0) and the computation performs as follows: 

/«L,./|/)|0)|0)H/)(^)|o)H/)( l">l°';^"l°' ) 

Then, by applying the CNOTqi; that is, the CNOT gate, defined by equations Ullfl . with 
the first qubit as the control one, we find: 



/ (g) CNOTqi 



f) 



0) I 0)+ I 1) I 0) 
V2 



f) 



0) I 0)+ I 1) I 1) 
V2 



f)\Poo), (19) 



where | /3oo) is the Bell state defined in expression Similarly, the Alice's circuit operates 
on the state given by expression 119() . The result has the general form: 



After -Alice — ^ Oij^fc | ijfc), 
ij,fc6{0,l} 



(20) 



Then, Alice measures the first and second qubits. Thus, the state just after the mea- 
surement will be one of the: 



I ^)o = aooo I 000) + aooi | 001); | ?A)2 = aioo | 100) + ami | 101); (21) 
I V)i =aoio I 010) +aoii I Oil); | V)3 = aiio | HO) + am | 111). 

Each result will be processed by Bob's part. If we trace each measurement result we can 
confirm that the initial state of the second qubit was delivered to the zeroth qubit, which 
belongs to Bob. For instance, let us suppose that the result was | ■0)0 given above. Then, 
Bob's circuit will outputs the following state j31| : 

(|0)|0)-|l)|0))(p|0)+g|l))^(|0)|0)-|l)|0))|/), 

which is a desired one, once the state of the second qubit was delivered to the zeroth; that 
is, it was teleported. The final state of the second qubit is different from the original one, 
which is accordance with the non-cloning theorem of quantum mechanics |19j . 

Once we have a quantum circuit (Figure 01 that performs the required computation, 
an interesting question arises: Given that quantum circuit how to find another one which 
performs the same computation but has less elementary gates? This optimization problem 
was addressed in through genetic algorithms (GAs). 

In fact, in |31|. a circuit for quantum teleportation is encoded by a chromosome that is 
a string of integers chosen from the set {0, 1, 2, 3} . Each gene is interpreted with a codon, 
i.e., a three-letter unit. The first letter indicates a kind of gate, whereas the second and the 
third letters indicate the qubits that the gate will operate. For instance, let us consider the 
following string: 

112 • 231 • 001 • 331 • 132 • 012 • 122 • 302 • 203 • 220 • 020 • 001 (22) 

The first codon whose first letter is 3 is interpreted as the partition between EPR-pair 
generation and Alice's part. The second codon whose first letter is 3 corresponds to Alice's 
measurement. The codification used in |3T] is defined by Tables |21 and |3| 








1 


2 3 







1 


2 3 





CNOTqi 


CNOTio 








CNOT12 


CNOT21 








CNOToi 


CNOTio 


1 





CNOT12 


CNOT21 


1 





CNOToi 


CNOTio 


2 





CNOT12 


CNOT21 


2 









3 









3 


1 


Lo 


Li 





1 


Li 


L2 





1 


Lo 


Li 


1 


1 


Li 


L2 


1 


1 


Lo 


Li 


2 


1 


Li 


L2 


2 


1 






3 


1 






3 


2 


Ro 


Ri 





2 


Ri 


R2 





2 


Rq 


Ri 


1 


2 


Ri 


R2 


1 


2 


Ro 


Ri 


2 


2 


Ri 


R2 


2 


2 






3 


2 






3 


3 






* 


3 






* 



Table 2: Codification for EPR generation (left). Alice's gates codification (right). 
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Table 3: Codons for Bob's circuit generation. 

This circuit representation is closer that one proposed in (described above). However, 
differently from that reference, in this representation the same integer may have different 
meaning in different parts of the chromosome. For example, the letter means CNOToi 
during the EPR-pair generation while it means CNOT12 for Alice's part. 

Once established the circuit representation, genetic operators based on mutations and 
crossover shall be specified. Mutations are implemented by properly change the alleles of a 
given chromosome. A two-point crossover is implemented by randomly choosing two parent 
chromosomes and exchange their alleles. 

Finally, the following steps are executed: (1) Decode each chromosome in a circuit and 
its gates; (2) Apply these transformation on the initial state; (3) If the circuit outputs a 
final state similar to the desired one, its fitness is enlarged. (4) Apply the genetic operators 
to generate the next population; (5) Go to step (1). 

Each individual (circuit) is evaluated as follows: 

(1) Make three random numbers a, /?, 7 G [0, 2n] . 

(2) Prepare three initial states (p, 9) given by: (e*'' cos a, e*'' sina) , 
(e*'*' cos /3, e*" sin (3) , (e*" cos 7, e''^ sin 7) 



(3) Derive the circuit by decoding the string using Tables |21 and |31 

(4) Use the state (p | 0) + g | 1)) ® | 00) as the input one. 

(5) Evaluate the fitness. 

(6) Change random numbers every 50 generations. 

After the measurement, we have to trace all the branches, corresponding to the possible 
outcomes given by equations H21(l . The desired final state, at the end of Bob's circuit has a 
general form (a | 00) + 6 | 01) + c | 10) + d | ll))«)(p | 0) +9 | 1)) = {ap,aq,bp,bq,cp,cq,dp,dq) . 
Thus, the ratio aij aij^x — pj q\ i = 0, 2, 4, 6. This gives a clue to find out an efhcient fitness 
function. 

Given a circuit, we observe that each one of the initial states produces 4 final states (one 
for each possible measurement's result, given by expressions (|21|l ). Thus, for the three initial 
states, we will have 12 possible final states. If we write the final states as (oj^o; ■■■7 0^,7), 
j — 0, 2, .., 11, we can express the gap between a final state and the desired one as: 

where n is the numbers of pairs such that (a^^i, flj i+i) 7^ (0, 0) and the summation is taken 
over such pairs. In the case that the final state is then error j is set to 100. The fitness 
function / is defined as: 



1 + 10 ^ error j 

If / is 1, that is, if the circuit is correct, the bonus of 1/ [number — of — gates) is added 
to / so as to apply a selection pressure based upon the circuit size (fitness is enlarged). 

In |2l] authors used the roulette- wheel selection and two-point crossover with probability 
0.7. Differently from the work found in for circuit design fsection f3.2ll . mutation is con- 
sidered more significant in this case. The mutation probability is 1/ [chromosome — length); 
that is, the algorithm is biased in the preservation of smaller chromosomes against the larger 
ones. The population size and the maximum number of generations was 5,000 and 1,000, 
respectively. All individuals are replaced every generation (there is no elitism) . The simpler 
circuit so obtained is pictured on Figure 
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Figure 5: Optimized circuit for the quantum teleportation problem. 

We can check that the circuit on Figure El is encoded by the string given by H22|l. This 
circuit has 8 gates while Brassard's had eleven (FigureQ)), which demonstrate the capabilities 
of the GA procedure to evolve an initial circuit towards a simpler one. 

We are in charge with an implementation of this algorithm, but using different strategies. 
We shall return to this point in sectional 



4 Quantum Evolutionary Computation 



4.1 Quantum-Inspired Genetic Algorithms 



In this section we present the second part of our review: the analysis of genetic algorithms 
based on quantum computing concepts. This is an important step towards the implemen- 
tation of genetic algorithms in a quantum hardware. 

We start with the Quantum-Inspired Genetic Algorithm (QIGA) proposed in JJ^. The 
QIGA is characterized by principles of quantum computing including qubits and probability 
amplitude. It uses a qubit representation instead of the usual binary, numeric, or symbolic 
representations ^jjTQ. More specifically, QIGA uses a m-qubit representation, defined as: 
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(23) 



where each pair (ctio, Pn) ; i = 1, m, indicates a qubit. 

Now, we must explain how convergence can be obtained with the qubit representation. 
Let us consider the following scheme, which is proposed in (11) . 

For each m-qubit chromosome of the form ^2'6\i . a binary string {xi, X2, Xm} is defined, 
where each bit is selected using the corresponding qubit probability, |aio|^ or jaiij^. Observe 
that if lofiol^ or |aii|^ approaches to 1 or 0, the qubit chromosome converges to a single state 
and the diversity given by the superposition of states disappears gradually. 

An application dependent fitness function is used to evaluate the solution {xi, X2, Xm}- 
Another step is to design efficient evolutionary strategies. This would be accomplished 
through crossover and mutations but their implementations are not explained in |11) . Ob- 
viously, as usual, we can suppose a one-point crossover between parent chromosomes as well 
as unitary operators to change a randomly chosen qubit {aio,l3ii) of the expression (|23|l . 
However, as the QIGA has diversity caused by the qubit representation, the role of genetic 
operators is not clear. Also, it is stated in jTJ that, if the probabilities of mutation and 
crossover are high, the performance of the QIGA can be decreased notably. 

At the beginning of the algorithm, a population Q (t) — {ql, 9m} of m-qubit chro- 

mosomes is instantiated. Given a m-qubit chromosome in Q (t) we can find the corresponding 
binary string through the rule stated above. The so obtained binary string population will 
be denoted by P (t) . 

Besides, there is an update step which aims to increase the probability of some states. 
Henceforth, given a qubit (aio^an) of a m-qubit chromosome, it is updated by using the 
rotation gate U (Oi): 
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(24) 



where 9i is formed through the binary solutions P (t) and the best solution found (see next) 
Let us present a pseudo-code of the QIGA developed in (llj : 
Procedure QIGA 
begin 

t ^ 

Initialize Q (t) 

Make P (t) by observing Q (t) 
Evaluate P (t) 

Store the best solution b among P (t) 
whilelnot termination-condition) do 
begin 

t^t+l 

Make P (t) by observing Q (t - 1) 



Evaluate P {t) 

Update Q [t] using quantum gates U {t) 
Store the best solution b among P (t) 

end 

end 

The quantum gates U (t) are application dependent. This step aims to improve the 
convergence. After updating Q (t) , the best solution among P (t) is selected, and if the 
solution is fitter than the stored best solution, the stored solution is replaced by the new 
one. The binary solutions P (t) are discarded at the end of the loop. A parallel version of 
the QIGA is presented in [T^ . 

A significant point to be considered is the exploration of the tensor product to enlarge 
diversity. Despite of authors claim in m, the scheme proposed did not take advantage of 
such effects at all. We will analyze this point in sectional 



4.1.1 Experiments for QIGA 

The knapsack problem, which is a kind of combinatorial optimization problem |18|. is used 
in to investigate the performance of QIGA. The 0—1 knapsack problem is described as: 
given a set of m items and a knapsack with limited capacity C, select a subset of the items 
so as to maximize a profit function / (x) given by: 



(25) 



and subjected to: 



J2^^^^<G, (26) 
i=l 

where (xi, 12, a^m) S {Oil}™; Pi and Wi are the profit and the weight associated to the 
item i, respectively. 

When applying the QIGA to this problem, the length of a qubit chromosome is the same 
as the number of items. The i — th item can be selected for the knapsack with probability 
jaiol^j following the procedure given in section r4.1l Thus, from each 771-qubit chromosome a 
binary string of the length m is formed. The binary string Xj represents the j — th candidate 
solution to the problem. The i — th item is selected for the knapsack if and only if Xij = 1. 
To measure the efficiency of the QIGA, its performance was compared with that one of 
conventional genetic algorithms (CGAs). Three types of CGAs were considered: algorithms 
based on penalty functions, algorithms based on repair methods, and algorithms based on 
decoders [TT|[^. 

For the first group of algorithms, the profit function is: 

m 

/(x) = ^p,x, -Pen(x), (27) 
1=1 

where Pen (x) is a penalty function. Among the possibilities to define penalty functions the 
following ones were considered in |11| : 



Peni (x) = log2 + P w^x, - ^ j j ' (28) 
Pen2 (x) = p lj2 ^'-^^ - ^ ) ' (29) 
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where p = max {pi/wi ; i=l,...,m}. 

For repair methods, the profit is defined by: 



m 



fix) = ^PrX[ 



(31) 



where x' is a repaired vector of the original vector x. Original chromosomes are replaced 
with a 5% probability in the experiment. The two repair algorithms considered in Jl] differ 
only in selection procedure, which chooses an item for removal from the knapsack: 

i?epi (random repair): The selection procedure selects a random element from the knap- 
sack. 

i?ep2 (greedy repair): All items in the knapsack are sorted in the decreasing order of their 
profit to weight ratios. The selection procedure always chooses the last item for deletion. 

A possible decoder for the knapsack problem is based on an integer representation. Each 
chromosome is a vector of m integers; the i — th component of the vector is an integer in the 
range from 1 to m — i + 1. The ordinal representation references a list L of items; a vector 
is decoded by selecting appropriate item from the current list. The two algorithms for this 
class used in ^T] are: 

ZJeci (random decoding): The build procedure creates a list L of items such that the 
order of items on the list corresponds to the order of items in the input file which is random. 

Z?ec2 (greedy decoding): The build procedure creates a list L of items in the decreasing 
order of their profit to weight ratios. 

Besides, there were an experiment that implements a scheme using Pen2 and Repi. 

The QIGA proposed in 11 for this problem contains a repair algorithm. It can be 
described as follows: 

Procedure QIGA-Knapsack 



t ^ 

Initialize Q {t) 

Make P (t) by observing Q {t) 
repair P [t] 
Evaluate P (t) 

Store the best solution b among P (t) 

while{t < MAX.GEN) do 

begin 

Make P (t) by observing Q {t ~ 1) 
repair P (t) 
Evaluate P (t) 

Update Q (t) using quantum gates U (t) 
Store the best solution b among P (t) 

end 

end 

Procedure make(x) 

begin 

i <- 

while(i < m) do 
begin 



begin 



if random [0, 1] > |aio|^ 

then <— 1 
else Xi ^ Q 

end 

end 

Procedure repair(a;) 

begin 

knapsack-overfilled ^ false 

if J2iLl WrXi > C 

then knapsack-overfilled <— true 
while (knapsack-overfilled) do 
begin 

select an i — th item from the knapsack 

if 1 'Wi^i < C' 

then knapsack-overfilled false 

end 

while (not knapsack-overfilled) do 
begin 

select a j — th item from the knapsack 
Xj ^ 1 

if 1 "^i^i > ^ 

then knapsack-overfilled <— true 

end 
^ 

end 

The profit of a binary solution is evaluated by expression and it is used to find 
the best solution h among P (t). A m-qubit chromosome is updated by using the rotation 
gates, following expression 1241) . The angles 9.i are computed as follows. Let us suppose that 
we have a binary string x = {xi,X2, ■■■,Xm) such that / {x) > f (5), where / is defined by 
expression (|31l) . If a:^ = 1 and bi = 0, the idea is to set the value of 9i = sign (ctio • an) A9i 
such that the probability amplitude of | 1) is increased. Thus, we want that \a'^i\^ > \aii\^, 
where a[i is given by equation H24|l . So: 

\a'ii\^ = {aMsinOi + an cos 9,)"^ = {aio sm9i)'^ + {an cos 9if + 2aioa,i sinft; cos6'i, (32) 

where we have supposed that aiQ,aii are real ones for simplicity. Thus, to increase the 
desired probability amplitude as much as possible we should set sign (aio ■ an) — +1, —1, 0, 
according to aioan > 0, aiQan < 0, or ai^an = 0, respectively. The setting of A6'i is 
through experimentation. In the reported example, it was set to 0.025. Following such 
procedure, a lookup table for 9i can be performed (see details). 

The update procedure is given bellow; 

Procedure update((7) 

begin 

i ^ 

while [i < m) do 
begin 

i^i + 1 

determine 9i 

obtain q' = {a^Q,a^i) as: 



(a-Q, a-i)^ = U {9i) (a^o, omf 

end 

q' 

end 

The results obtained by the QIGA just presented, reported in fTTj, uses the following 
profits and weights: 

Wi — uniformly — random[l, 10), (33) 
pi — Wi + 5. 

The average knapsack capacity was used: 

C^-Y^w,. (34) 

1=1 

The data files were unsorted and the number of items were 100, 250 and 500. 

The population size of the eight conventional genetic algorithms was equal to 100. Prob- 
abilities of crossover and mutation were fixed: 0.65 and 0.05, respectively. The population 
size is 1, for the first series of experiments, and 10 for the second one. As a performance 
measure of the algorithm the best solution found within 500 generations over 25 runs is 
collected. Also, the elapsed time per one run is checked. 

For 100 items QIGAs yielded superior results as compared to all the CGAs. For 250 and 
500 items the QIGA with 10-size population outperforms all the classical ones |11| . 

4.2 Quantum Genetic Algorithms 

The work reported on section ETI shows that the application of quantum computing concepts 
to evolutionary programming is a promising research. The results presented points out that 
a quantum genetic algorithm (QGA) would outperform the classical ones. Besides, such 
implementation would take advantage of quantum parallelism as well as GAs parallelism. 
The obvious question is how to implement genetic algorithms in quantum computers? 

The reference is an effort to produce a QGA. Despite of the fact that there are 
several open points, it is the first effort in the direction of such algorithm. 

The QGA proposed in |25| uses two registers for each quantum individual; the 
first one stores an individual while the second one stores the individual's fitness. 
These two registers are referred as individual register and the fitness register, respec- 
tively. A population of N quantum individuals is stored through pairs of registers 
{individual — registeri, fitness — registeri), i = 1, .., N. 

At different times during the QGA the fitness register would store a single fitness value 
or a quantum superposition of fitness values. Identically for the individual register. 

Once a new population is generated, the fitness for each individual would be calculated 
and the result stored in the individual 's fitness register. 

The effect of the fitness measurement is a collapse given by expression © . This process 
reduces each quantum individual to a superposition of classical individuals with a common 
fitness. It is a key step in the QGA |25| . Then, crossover and mutation would be applied. 
The whole algorithm can be written as follows: 

Quantum Genetic Algorithm 

Generate a population of quantum individuals. 

Calculate the fitness of the individuals. 

Measure the fitness of each individual (collapse). 

■u;ft,iZe(termination-condition) do 

Selection based on the observed fitness. 



Crossover and Mutations are applied. 
Calculate tlie fitness of tlie individuals. 
Measure the fitness of each individual (collapse). 

end while 

According to the more significant advantage of QGA's will be an increase in the 
production of good building blocks {schemata llHj ) because, during the crossover, the 
building block is crossed with a superposition of many individuals instead of with only one 
in the classical GAs. 

One can also view the evolutionary process as a dynamic map in which populations tend 
to converge on fixed points in the population space. From this viewpoint the advantage of 
QGA is that the large effective size allows the population to sample from more basins of 
attraction. Thus, it is much more likely that the population will include members in the 
basins of attraction for the higher fitness solutions. 

Another advantage is the quantum computer's ability to generate true random numbers. 
By applicating Kolmogorov complexity analysis, it has been shown that the output of classi- 
cal implementations in genetic programming, which use a pseudo random number generator, 
are bounded above by the genetic programming itself, whereas with the benefit of a true 
random number generator there is no such bound |25l I24j . 

Despite of these promising features, fundamental points are not addressed in Firstly, 
it is not clear how to implement crossover in a quantum computers. Besides, how to per- 
form the fitness function calculation in quantum hardware? Even a much more funda- 
mental problem is that to explore the superposition of quantum individuals the correlation 
individual fitness must be kept during the whole computation. Entanglement seems 
to be the only possibility to accomplish this task. But, in this case, things must be for- 
mally described to avoids misunderstandings and wrong interpretations. We develop such 
mathematical formalism on section |S| 

5 Discussion and Perspectives 

In this section we analyze some issues concerning to the reviewed methods. Possible solutions 
and perspectives in this area are also discussed. 

In iBJ we show some challenges concerning the GA for learning linear operators (section 
\'^.^ . Other tests presented in "5" show that the number of generations seems to increase 
when space dimension gets higher. The increasing rate must be controlled if we change 
the population size properly. However, such procedure could be a serious limitation of the 
algorithm for large linear systems. 

The behavior for underconstrained problems; that is, when S = 
{(I Xi), I "^i)) ; P\Xi)=\^i), i = 1, K < n} , where n is the space dimension, is 
also analyzed in ^\ In this case, we had to increase the population size but the number of 
generations is smaller than that one for the constrained test {K — n). As we expect, there 
is a trade-off between the increase of solutions and the fact that we are less able to properly 
evolve the populations due to the lack of prior information. Moreover, the observed success 
is an advantage of the method, if compared with traditional ones. In this case, numerical 
approaches based on iterative methods in matrix theory (Gauss-Seidel, GMRES, etc) can 
not be applied without extra machinery because the solution is not unique 'H'. 

The comparison with Dan Ventura's learning method, given on section Td.l.ll shows that 
our algorithm overcomes the limitation of the later: we do not need that {| xo), I Xi)} is an 
orthonormal basis of the vector space. However, when using our GA method, we pay a price 
due to storage requirements and computational complexity. 

Dan Ventura's algorithm as well as numerical methods (see |H] and references therein), 
have a computational cost asymptotically limited by O (n * n^) while our GA method needs 



{Ngen * N * n^) float point operations. Besides, for traditional numerical methods and 
Ventura's algorithm, we observe a storage requirements of O (n^) against O [N * 71^) for our 
approach. Thus, the disadvantage of our method becomes clear. 

However, if compared with matrix methods, our algorithm is in general less sensitive to 
roudofF errors '10'. This is due to, unlikely numerical methods that try to follow a path 
linking the initial position to the optimum, our GA algorithm searches the solution through 
a set of candidates. 

To improve the convergence we need better evolutionary (crossover/mutation) strategies. 
The behavior pictured on Figure |21 of section 15. 1.11 is a typical one for every test we made 
It indicates that our evolutionary strategies are efficient to get closer the solution but 
not to complete the learning process. Further analysis should be made to improve these 
operators. 

When comparing the works for quantum circuit design, |2Hj and |H1| . we observe the 
following aspects: 

1) Gate representation: The gate structure of versus the codon used in PJ. Despite 
of some apparent difference between them, it is simple to check that they are equivalent, in 
the sense that, any gate can be represented with either the later or the former. 

2) Genetic Operators: Both implementations have used crossover and mutations. How- 
ever, the crossover implementation used in |23| operates on all the levels of an individuals 
structure (the gates, each category of qubit operands and each parameter type) while in '31j 
it affects only the gate and qubit levels. Mutations are basically equivalent because, if the 
alleles are randomly changed, like in |31] . we are randomly replacing a gate with a new one, 
like in |28j . and vice- versa. 

3) Range of Applications: Despite of the fact that the aim of [2] is circuit optimiza- 
tion, it can be straightforwardly adapted for circuit design. This can be accomplished by 
changing the mutation probability (the formula 1/ {chromosome — length) does not make 
sense in this case). We can follow [53] and set this probability to a small value (typically 
0.001). Besides, the fitness function remains case dependent and we stop using the bonus 

1 / (number — of — gates) to bias the solution to smaller circuits (if we do not know any 
prior correct circuit, there is no sense for prefering smaller circuits over bigger ones during 
evolution). Besides, some elitism may be introduced. Now, we are analyzing such modifica- 
tions. 

Moreover, a more fundamental question about GAs for circuit design follows from the 
next comments. The algorithms |23l 13 1| basically evolve an initial population of individuals 
towards a desired circuit. Evolution can be regarded as the exploration of search spaces by 
populations. Thus, an interesting question in this case is what about the structure of such 
spaces for circuit designing/optimization? 

For instance, we must observe that without the identity operator the circuit size (number 
of gates) will be a variable. Henceforth, a search space of 12-codons strings (like expression 
H22|ll. would be transformed in another one with just 8-codons strings at the end of the 
optimization process. This can be seen as an evolutionary process called innovation 2(1 . 

Following |26| we do need a mathematical representation in which the kind and number 
of codons follow from the dynamics of the model. In [2H the concept of configuration spaces 
is proposed as one of such approach. Obviously, the Identity operator is a simple trick to 
address this problem if we know in advance the maximum circuit size. However, the concept 
of configuration spaces might open possibilities to analyze the structure of the circuit space. 
That is way we are going to consider this mathematical framework in our research. 

A configuration space is a set of objects (circuits, for example) as well as a topological 
structure on this set which describes how these objects can be transformed into each other 
by an operator |26' . 

Symmetries of the configuration space induced by evolutionary mechanisms (mutation 
and crossover, for instance) are fundamental elements in this framework. They define the 



dimensionality of the space in which evohition occurs. Hence, any evolutionary process that 
affects the symmetries of the configuration space may change its dimensionality (the number 
of non- identity gates, in our case). The configuration space formalism includes beautiful 
mathematical results in finite Abelian groups |26j . We wish to analyze the algorithms 
proposed in |23[ 1311 through this formalism, in order to derive more efficient evolutionary 
strategies. 

When considering the QIGA presented on section l^lTl some explanations must be offered 
about the diversity that can be achieved by the m-qubit representation .11 . 

Thus, let us take the make(x) procedure. For simplicity, we are going to consider a 
3-qubit chromosome. Let ri,r2,r3 be the random numbers generated during the while loop 
execution in make(a;). We could have ri > |q!io|^,''2 > |a2ol^;''3 < |q^3oI^- Thus, the 
generated string would be (1,1,0). 

Now, consider the tensor product 

I Vi)® I '02)«) I V'a) = X! aina2»2a3J3 hi>® I «2)«) I is)- (35) 

il,i2,'i36{0,l} 

Thus, the qubit chromosome will be represented as a superposition of the states | | 
*2)<8' I is), ?i,*2j*3 G {0, 1}, and so it carries information about all of them at the same 
time. Such observation points out the fact that the qubit representation has a better charac- 
teristic of diversity than classical approaches, since it can represent superposition of states. 
In classical representations we will need at least 2^^ chromosomes to keep the information 
carried by expression (|35|l while only one 3-qubit chromosome is enough. 

However, the probability amplitude of the state | 1)(8' | l)^ | 0) may not be the largest 
one. Henceforth, it does not seems that the binary string generation rule proposed in |11| 
does explore such diversity in general. 

However, if we take another generation rule, say: if jaal^ > l^iol^ then Xi <— 1, else 
Xi <— 0, thus we can be sure that the generated binary string is an index to the larger 
amplitude probability of state 1)35(1 . Experiments must be performed in order to show the 
efficiency of such rule. 

The QGA presented in !25! exploits the quantum effects of superposition and entan- 
glement. However, the lack of a more formal explanation raises some questions. How to 
implement crossover in quantum computers? How to compute the fitness function? What 
about a mathematical definition of a quantum individual? These are examples of such 
questions. 

Now, we address some of these points in order to be closer to answer the question: what 
GAs will look like as an implementation on quantum hardware? 

The starting point of our development comes from the known problem of finding the 
period r of a periodic function / : Zn Z, where Zn denotes the additive group of 
integers modulo N. 

In this case, the quantum solution provided by Shor |19j uses a hardware with two 
registers in the following entangled state: 

N-l 
^ x—O 

Thus, according to the expression ®, by measuring the second register, yielding, say, a 
value uo, the first register 's state will collapse to an uniform superposition of all those \ x)'s 
such that / (x) = yo; that is: 

1 ^"^ 

l*>a/ter = ^5Z \ ^0 + kr) , (37) 



where xq is such a x and N = Kr. When using the state | ^'), given by expression (|35|) . the 
desired effect is to keep the correspondence between each integer x with its corresponding 
value / {x) . 

Now, let us return to the QGA of section W?^ and present a physical description of it. 
The quantum individual could be mathematically represented by a state given by expression 
H36|l ■ where | x) represents an individual and / (x) its fitness. Thus, we keep the idea of 
representing a quantum individual through two registers which was used in So, if 

we have M quantum individuals in each generation we need M register pairs {individual 
register, fitness register). 

In our formulation, each register is a closed quantum system. Thus all of them can be 
initialized with the state given by expression H3t)|l . Then, unitary operators W will be applied 
in order to complete the generation of the initial population. Henceforth, the initialization 
could encompass the following steps: 

1) For each register i, generate the state: 

^ a:— 

2) Apply unitary operators W (rotations, for example) and [//, the known black box 
which performs the operation Uj \ a)® \ 0) =| a)® \ f (a)) to complete the initial 
population: 

I VI/), ^ UfW I = Uf (w (^^) ® I 0), 




N-1 

= J2^f ("^ I ^^^"^ I = 

x=0 

N-l 

= 51 I ^)^® I /(^))''' ^^'^,-,M. (38) 

We must highlight that all the above operations are unitary ones, consequently, can be 
performed in quantum computers [201 ■ Besides, it is important to observe that the fitness 
is stored in the second register after the generation of the population. Now, by measuring 
the fitness, each individual undergoes collapse, according to the expression (|37|l : 

where | k)i is such that the observed fitness for the i — th register is / {k) = j/q. 

When entering the main loop, the observed fitness is used to select the best individuals. 
Then, genetic operators must be applied. 

Mutations can be implemented through the following steps. 

1) Apply Uy^ over the measurement result: 

Uyl\^yj^er^lj2 |%®|0)„ (40) 
k=0 

2) Unitary operators P (small rotations, for example) are applied to the above result: 



P {uj' I *>r^"'^) - E ^ (1^) ® I 0). = 

2:=0 

where we expanded the resuh in the computational basis. 

3) Finally, apply Uf to recover the diversity that was lost during the measurement: 

N-l 

UfPUy' I vi/)^/*- = E I ^)^® I / (42) 

x=0 

The development given above allows to discuss some points. Firstly, we observe that 
if we take a superposition of individuals in the first register and the corresponding fitness 
superposition in the second one, as claimed in I25[ . we will have: 

(N-l \ /N-l \ Af-lAf-1 

Y.a,\x)UiY,b, \f{x)) = E E ^-^y I I /(y))- 
x=0 ) \x=a ) x=a y=a 

Thus, we are not able to keep the correlation individual <-> fitness. For instance, after 
a measurement that gives a zq value, the system state would be: 

N-l 
x=0 

where / (yo) = zq (observe that in general x yo in this expression). So, such proposal does 
not seems to be efficient at all. 

According to 25 , the major advantage for a QGA is the increased diversity of a quantum 
population due to superposition, which we have precisely defined through expression (|38[l . 
This effective size decreases during the measurement of the fitness, when the superposition is 
reduced to only individuals with the observed fitness, according to expression 1)39(1 . However, 
it would be increased during the crossover and mutation applications. Besides, by increasing 
diversity it is much more likely that the population will include members in the basins of 
attraction for the higher fitness solutions. Thus, an improved convergence rate must be 
expected. Besides, classical individuals with high fitness can be relatively incompatible; 
that is that any crossover between them is unlikely to produce a very fit offspring. However, 
in the QGA, these individuals can co-exist in a superposition. 

Despite of the mathematical development given above, two fundamental points remain. 
Firstly, we can not suppose that the number of elements of the search space is the same of 
the number of states of the computational basis (N, in the above presentation). If so, the 
solution would be to find the maximum value of / (x), which is just an optimization problem 
that can be addressed by quantum optimization algorithms |22|- Besides, the search space 
size is in general too large that makes some assumption unreasonable. 

Secondly, the crossover needs special considerations not only because combination of 
states in Hilbert spaces is limited by the constraint of unitary operations but also because 
each register pair is a closed quantum system. Thus, we need some kind of quantum com- 
munication channel to combine states. This question should be addressed in the context of 
state-of-the-art quantum computers architecture (see |2(J| and references therein). 



6 Conclusions 



In this paper we survey the main works in quantum evolutionary programming and in the 
apphcations of GAs to address some problems in quantum computation. Besides, we offer 
new perspectives in the area which are part of our current research in this field. Among them, 
we believe that the analysis of the algorithms proposed in [531 1^ through the configuration 
space formalism and a QGA implementation are the most exciting ones. 

The concept of configuration spaces might open possibilities to analyze the structure of 
the circuit space in order to derive more efficient evolutionary strategies. 

On the other hand, a QGA implementation could take advantage of both the quantum 
computing and GAs parallelism. We analyze the work summarized on section Pl.2l and give a 
formal explanation of its main elements. However, quantum crossover and efficient strategies 
for search space exploration remains challenges in this field. 
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